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ABSTRACT 

We measure the star formation rate (SFR) as a function of redshift z up to z ss 4.5, based on B, I 
and (I+B) selected galaxy catalogues from the FORS Deep Field (FDF) and the K-selected catalogue 
from the GOODS-South field. Distances are computed from spectroscopically calibrated photometric 
redshifts accurate to l^z/{zspec + 1) < 0.03 for the FDF and < 0.056 for the GOODS-South field. 
The SFRs are derived from the luminosities at 1500 A. We find that the total SFR estimates derived 
from B, I and I+B catalogues agree very well (< 0.1 dex) while the SFR from the K catalogue is 
lower by « 0.2 dex. We show that the latter is solely due to the lower star-forming activity of K- 
selected intermediate and low luminosity [L < L^) galaxies. The SFR of bright [L > L*) galaxies 
is independent of the selection band, i.e. the same for B, I, (I-|-B), and K-selected galaxy samples. 
At all redshifts, luminous galaxies {L > L*) contribute only ~ i to the total SFR. There is no 
evidence for significant cosmic variance between the SFRs in the FDF and GOODs-South field, < 0.1 
dex, consistent with theoretical expectations. The SFRs derived here are in excellent agreement with 
previous measurements provided we assume the same faint-end slope of the luminosity function as 
previous works {a ^ —1.6). However, our deep FDF data indicate a shallower slope of a = —1.07, 
implying a SFR lower by w 0.3 dex. We find the SFR to be roughly constant up to z « 4 and then 
to decline slowly beyond, if dust extinctions are assumed to be constant with redshift. 
Subject headings: galaxies: high-redshift — galaxies: formation — galaxies: evolution 



1. INTRODUCTION 

The determination of the star formation rate (SFR) 
history of the universe is one of the most interesting 
results extracted from the deep photometric and spec- 
troscopic surveys of the last decade. A large num- 
ber of measurements have been collected, at low (the 
Cana da-France redshift survey at z < 1, iLillv et aP 
11596'), and high redshift from the Hubble Deep Field 
North (Madau ct al. 1996) , the large samp les of U and 
B drop-out galaxies IjSteidel et alJ Il999ji . up to the 
most recent determinations based on I-dropouts at red- 
shift w 6 f rom the GOODS, UDF and UDF-Parallel 
ACS field s lIGiavalisco et al.ll2004al IBunker et al1l200l 
IBouwcus e t al.ll2004j) . These studies show that the SFR 
(uncorrected for dust) increases from z = to z = 
1, stays approximately constant in the redshift range 
1 — 4, and starts to decline at larger redshifts. In all 
the cases quoted above the determination is based on 
the estimate of the total UV galaxy luminosity den- 
sity, that for a given IMF is prop ortional to the instan- 
taneous SFR ijMadau et all Il99a . 1998 , ) As discussed 
by many authors (e.g. iHopkins et al I l20?Ul) this ap- 

^ Based on observations collected at the European Southern Ob- 
servatory, Chile (ESO programmes 63.O-0005, 64.O-0149, 64. 0- 
0158, 64.O-0229, 64.P-0150, 65.O-0048, 65.O-0049, 66.A-0547, 
68.A-0013, and 69.A-0014) 



proach is affected by the uncertainties of dust correc- 
tion, but roughly agrees with other estimators at low 
to intermediate redshifts (z < 1). Theoretical mod- 
els of galaxy formation and evolutio n can be tested 
against the measured SFR history ijSomerville et alJ 
iOOl; Hcrnquist & Springcl 2003). 

So far, all determinations of the SFR history have suf- 
fered from some major limitations. High redshift samples 
have been small in number due to the limited field of view 
of deep pencil-beam surveys, resulting in large Poisso- 
nian fluctuations and large field-to-field variations (cos- 
mic variance). The faint-end of the luminosity function 
(LF) is thusfar only poorly constrained at high redshifts, 
implying large completeness correction factors. Finally, 
the technique used to generate the high-redshift galaxy 
catalogues (drop-out selection, optical magnitude limited 
survey) might have introduced biases by selecting only 
specific types of gala xies and possibly missing relevant 
fractions of UV hght l|Ilbert et alJl2004|) . 

Here we try to minimize these uncertainties and de- 
termine the SFR history of the universe with improved 
accuracy up to z w 4.5. Our sample of high redshift 
galaxies is based on two d eep fields, the (I an d B selected) 
FORS Deep Field (FDF, IHeidt et a,l.ll2f)f)l. and the (K - 
selected) GOODS-South field lIGiavahsco et al.ll2004bt) . 
Both cover a relatively large sky area, reducing the prob- 
lem of cosmic variance. Both are deep enough to allow 
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the detection of several x 10"^ galaxies, thus minimizing 
the effect of shot noise. 

Accurate photometric redshifts {/^z / {zspec + 1) < 0.03 
for the FDF and < 0.056 for the GOODS-South field) 
with only w 1 % catastrophic failures allow us to mea- 
sure the UV luminosity function down to fainter limits 
than spectroscopic samples. A detailed comparison of the 
UV luminosity functions of the FDF with the LF derived 
in large surveys was presented in IHabasch et al. ( 2004) 
and shows good agreement in the overlapping magnitude 
range at all redshifts. Finally and most importantly, the 
two fields provide us with B-band, I-band and K-band 
selected catalogues, making it possible to assess the de- 
pendence of the SFR on the detection band and galaxy 
colors and the associated selection biases. 

The Letter is organized as follows. In [J21 we dis- 
cuss the photometry and the photometric redshifts of 
the two fields, in i]31we present our results on the SFR 
history, and in 2] we draw our conclusions. Through- 
out the paper we use AB magnitudes and adopt a 
concordance cosmology with VLm = 0.3, f^A = 0.7, and 
Ho = 70kms-iMpc"\ 

2. DATA SETS 

The present results are based on photo i netric 
catalogues derived for the FDF ijHeidt et al l |2()0:it 
Eabasch ct al. 2004. UBgRI,8 34 nm,zJ Ks bands) and 
the GOODS-South fields (.Salvato et alJ[2004 UBVRI- 
JHKs bands). The two fields cover approximately the 
same area (39.81 arcmin^ for FDF and 50 arcmin^ for 
GOODS); the FDF reaches effective absolute magnitude 
limits « 1 mag deeper than GOODS-South (see below). 
We use the I-b and and B-b ar id sel ected FDF catalogues 
as de rived in iHeidt etlHI l|2flfl3tl and iGabasch et al.l 
l)2004jl . The B and I selected catalogues list 5488 and 
5557 bona-fide galaxies (having excluded the one known 
bright quasar in the field) down to Bum = 27.6 and 
Iiim = 26.8, respectively. The I-l-B catalogue obtained 
by merging these two contains 6756 entries. Photomet- 
ric redshifts for the FDF galaxies have an accuracy of 
I^z/i z.^ner + 1) < 0.03 wlt h Only ~ 1% catastrophic out- 
liers (jGabasch et al.ll2004() . 

Our K-band selected catalogue for the GOODS-South 
field is based on the 8 2.5 x 2.5 arcmin^ J, H, Ks VLT- 
ISAAC images publically available, taken with seeing 
in the range 0.4" — 0.5". The U and I images are 
from GOODS/EIS public survey, while B V R are taken 
from the Garc hing-Bonn Deep Surv ey. Data reduction 
is desc ribed in lArnouts et alJ l)2001|) and Schirmc r et aP 
l|2003|l . respectively. The data for the GOODS field 
were analyzed in a very similar way to the data of the 
FDF. The objects were detected in the K-band images 
closely following th e procedure used for the FDF I and 
B band detection llHeidt et alJl2003D . using both SEx- 
tractor ('Bertin & Arnouts^'19961) and the YODA pack- 
age (^rorv 2003). A detailed desc r iption of the proce- 
dure can be found in lSalvato et al.l 1)20041) . We detected 
3367 objects in K-band for which we derived magni- 
tudes (fixed aperture and total) in all bands. Number 
counts match the literature values down to K k, 25.4, 
which is the completeness limit of the catal ogue, in agree- 
ment with the number obtained following Snigula ct alj 
(|2002|) . Note that much deeper AC S based catalogues 
are available ijGiavalisco et aljf2004bj) . but as we are fo- 



cusing on the K-selection they are not relevant in this 
context. We comp uted photometric redshifts following 
iBender et aP pOOl 2004, in preparation) and using the 
same SED template spectra as for the FDF. The compar- 
ison with the spectros copic redshifts of the VIMOS team 
l|Le Fevre et al.ll200l) and the F0RS2 spectra released 
at http://www.eso.org/science/goods/, shows that the 
photometric redshifts have an accuracy l^z / [zspec-'r 1) < 
0.056. Similar results are obtain ed when comparing to 
the COMBO-17 (|Wolf et al.""200 4D data This is nearly a 
factor of 2 better than iMob asher et al.l l|2b04^) obtained 
using ground-based plus HST/ ACS data. Sta rs are iden- 
tified and exc hided as inlGabasch et aP l|2004D . as well as 
known AGN l|Szokolv et alJ2003|) . leaving 3297 bona-fide 
galaxies used in the further analysis. 

Fig. n shows the distribution of galaxies (slightly 
smoothed with a Gaussian kernel) in the rest frame 
1500A absolute magnitude M1500 vs. redshift plane, 
computed by integrating the best fitting SED over the 
band definition (1500 ± lOOA). The contours agree re- 
markably well at the bright-end showing that the number 
density of bright galaxies does not significantly depend 
on the wavelength at which they were selected (for the B- 
band, this is of course only true up to z ~ 3). For better 
comparability of the FDF and GOODS-South samples 
at the faint-end, we chose a consistent magnitude cut-off 
for all samples in Fig. Q This magnitude cut-off corre- 
sponds to the completeness limit of the GOODS-South 
sample and is about one magnitude brighter than the 
completeness limits of the FDF B and I samples. For the 
redshift bins defined by the limits 0.45, 0.81, 1.21, 1.61, 
2.43, 3.01, 4.01, 5.01, the cut-offs in Afisoo are at -15, 
-16, -17, -18, -19, -20, -20. 

The M1500 LFs and the related parameters M*, 
and a of the B, I, and (I-fB) selected FDF samples 
are almost identical. We derived consistent values for 
a (—1.07 ±0.04) for all th ree samples considere d here, 
similar to that described in lGabasch "eTall ll200l . Con- 
sistent faint-end slopes (a = —1.01 ±0.08) were obtained 
using a brighter subset of the data set (i.e. 1 mag brighter 
than the 50% completeness limit). Objects that were de- 
tected in only one band and not in both are all faint 
and do not contribute signi ficantly to SFR determined 
from the integral over the LF. IGabasch et aT] l)2004l) show 
that the steeper slope of other surveys is largely due to 
shallower limiting magnitudes. This is supported by an 
analysis of z ^ 6 dropouts fr om GOODS an d the Hubble 
Ultra Deep Parallel Fields l)Bouwens et al.ll2(T04) where 
an a of -1.15 was derived. Compared to LF parame- 
ters of the optically selected samples, the M1500 LF of 
the K-selected sample has slightly brighter values of A/, , 
significantly lower values in and, within the large er- 
rors, a similar faint-end slope a. Since the slightly shal- 
lower K-selected sample does not allow us to constrain 
the faint-end slope to the same level as our FDF sample 
(but is consistent with the faint end slope a ~ —1.07 
determined for that field), we adopt this value for our 
K-selected sample. 

We examine the consequences of these findings for the 
SFR in the next section. 

3. THE STAR FORMATION RATE 

We compute the SFR for all 4 catalogues from the 
total luminosity densities ^1500 in the 1500 A band. 



The SFR history in the FDF and GOODS-South fields 



3 



First, we derive Z1500 at a given redshift by summing 
the completeness correc ted (using a ViVmax correction, 
see'Gabasch e t al.ll200^ LFs up to the 1500 A absolute 
magnitude limits. Second, we apply a further correc- 
tion (to zero galaxy luminosity) ZGL, to take into ac- 
count the missing contribution to the luminosity density 
of the fainter galaxies. To this end we use the best-fitting 
Schechter function. For the FDF catalogues the ZGL 
corrections are only 2-20% in size. The small ZGL cor- 
rection employed here owes itself to the faint magnitude 
limits probed by our deep FDF data set and the relatively 
flat slopes (a « —1.07) of the Schechter function. Due 
to the brighter magnitude limit, the ZGL corrections for 
the GOODS catal ogue can be as hig h as 50%. Note that 
if we follow i.e. ISteidel et al.l il999|) who find a = —1. 6 
(excluded at 2ct with our fits, see iGabasch et ani2004D . 
we would get much larger ZGL corrections for the same 
M*, (see the dotted line in Fig. Inland the discussion 
below). 

Finally, following Madau et al. (1998) we derive the 
SFR by scaling the UV luminosity densities: SFR1500 — 
1.25 X 10^^^ X /i5oo in units of MQyr~^ M-pc'^ , where the 
constant is computed for a Salpeter IMF. The resulting 
values of SFR1500 are shown in Fig. [3 as a function of 
redshift. Errors are computed from Monte Carlo simu- 
lations that take into account the probability distribu- 
tions of photometric redshifts and the Poissonian error 
(IGabasch et al.l l2004|) . Following lAdelbereer fc Steidell 
l^200^)^ . we assume that dust extinction does not evolve 
with redshift and is about a factor of ~ 5 — 9 in the rest- 
frame UV. A more detailed discussion of the role of dust 
will be given in a future paper, like an analysis based 
on the SFR derived at 2800 A. Thanks to the large area 
covered and the faint limiting magnitudes probed, our 
determination of the SFR is the most precise to date, 
with statistical errors less than 0.1 dex for the single 
redshift bins spanning the range 0.5 < z < 5. 

The considerations of f|21translate in the following con- 
clusions about the SFR. Out to redshift z « 3 the SFRs 
derived from the I and B selected FDF, or the merged 
I-l-B catalogue, are identical within the errors (< 0.1 dex; 
see plot at the bottom left of Fig. [21 . At larger redshifts 
the B-selected SFRs underestimate the true values, since 
B drop-outs are not taken into account. The strong evo- 
lution in both the M, and parameters of the Sc hechter 
LF m easured as a function of redshift bv .Gabasch et al.l 
l|2004ll results in a nearly constant SFR, because the 
strong brightening of M, is compensated by the dramatic 
decrease of with z. Comparing the two lower panels 
of Fig. |2] shows that luminous galaxies (L > L*) con- 
tribute only a third of the total SFR at all observed z, 
independent of the selection band. 

The K-selected SFRs are similar in shape, but system- 
atically lower by « 0.2 dex at z > 1. This result holds 
independently of our completeness correction. If we con- 
sider only the contributions to the SFR down to the lim- 
iting magnitude set by the K-band, we find the same 0.2 
dex difference for 1 < z < 3, and 0.15 dex at z > 3. 
Fig. n shows that this result originates from the lower 
density of M1500 > — 19 galaxies in the K-selected cata- 
logue, as intermediate and low luminosity blue galaxies 
contributing to the SFR budget are more easily detected 
in the bluer bands than in K. In fact, the contributions 



to the SFR coming from the galaxies brighter than Ll 
are identical within the errors for the I and K selected 
catalogues (see Fig. |21 bottom-right panel). Therefore, 
cosmic variance does not play a role, as we also verified 
by comparing the B-band number counts in the 2 fields. 
They agree withi n 0.1 dex, which is the expected vari- 
ation derived by iSomerville et al.l l|2004ll scaled to the 
area of the GOODS -South field. On the other hand, 
IGabasch et~aLl l|2004D show that the I-band FDF cata- 
logue might be missing only about 10 % of the galaxies 
that would be detected in a deep K-band select ed sur- 
vey w ith magnitude limit Kab ~ 26 (like in Labb e et alJ 
iSl- The missing galaxies would be faint and likely not 
contributing significantly to the SFR provided their dust 
extinction is not exceedingly large. Independent of the 
selection band the SFR declines beyond z ~ 4.5. Our 
results confirm the conjecture of Kashikawa et al. (200^1) 
that the K-selected UV LFs match the optically selected 
LFs at high luminosities. 

The comparison with the literature shows that our re- 
sults are ^ 0.3 dex lower, independent of the selection 
band. This difference stems from the large com pleteness 
corrections apphed by, e.g.. ISteidel et al.l l|1999D . derived 
from the steep slopes fitted to the LF (see Our re- 
sults scale to the literature values if similar slopes are 
used for the same and 0, . This is shown by the dot- 
ted lines of Fig. El where we have assumed a slope of 
— 1.6 for our data set while keeping and the same 
as in our fit. 

The overall agreement between the SFRs derived over a 
wide wavelength range (within 0.2 dex), from the optical 
B and I to the NIR K, sampling at z ^ A the rest-frame 
UV and B, shows that we are approaching (in the opti- 
cal) the complete census of the galaxies contributing to 
the stellar production of the universe up to this redshift. 
Therefore, we can expect possible biases indu ced by miss- 
ing st ellar energy distributions with redshift l|llbert et alJ 
120041) to be small, when deep enough optical or NIR 
catalogues are available. However, we might still not 
take into account the possible contribution to the SFR 
coming fr om faint, high ly dust -absorbed r ed star-forming 
galaxies (iHughes et anil998t iGenzel et all 1200 iD which 
are likely missing from optically or near-infrared selected 
samples. Nevertheless, i t is encouraging to find that re- 
cent Spitzer results (e.g. lEgami et al.l200^ indicate that 
the majority of the star formation has already been ac- 
counted for using the dust-corrected SFR derived from 
optical studies. 

4. CONCLUSIONS 

We have measured the SFR of the universe out to 
z « 4.5 with unprecedented accuracy from the FORS 
Deep Field and the GOODS-South Field (90 arcmin^ in 
total). Our main conclusions are: 

• The cosmic variance in the SFR history between the 
FDF and GOODS-South field is negligibly smafl. The 
difference between these fields is < 0.1 dex, consistent 
with theoretical expectations. 

• The SFR of galaxies brighter than is the same (< 0.1 
dex) in B, I, (I-l-B) and K selected catalogues. This indi- 
cates that present optical and NIR surveys are unlikely 
to have missed a substantial fraction population of mas- 
sive star forming objects (with the possible exception of 
heavily dust-enshrouded starbursts). 
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• The total SFR integrated over all galaxy luminosities 
is the same in the B, I, and (I+B) selected catalogues 
and is lower in the K-selected catalogue by 0.2 dex. This 
difference originates at luminosities lower than which 
implies that K-selected surveys miss a significant fraction 
of star- forming lower-luminosity galaxies. 

• At all redshifts, luminous galaxies (L > L*) contribute 
only i to the total SFR, i.e. the integrated SFR of 
L < L^, galaxies is a factor of ~ 2 higher than the one of 
L > L:t galaxies. 

• Our fits to the FDF luminosity functions suggest a flat 
faint-end slope of a = —1.07 ±0.04 in contrast to the as- 
sumed slope of a —1.6 in the literature. This implies 
that past determinations have overestimated the SFR by 
a factor 2. 

• The SFR is approximately constant over the redshift 
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Fig. 1. — The distribution of galaxies in the rest-frame, 1500 A absolute magnitude vs. redshift space, slightly smoothed with a Gaussian 
kernel. The red colors refer to the K-selected galaxies of the GOODS-South field, the blue and the green colors to B and I selected galaxies 
of the FDF. The lowest contour corresponds to 0.75 galaxies/arcmin^/mag per unit redshift bin; the others give the 2.5, 3.75, 6.25, 8.75, 
11.25 and 13.75 galaxies/arcmin'^/mag per unit redshift bin density levels. For a better comparison of the FDF and GOODS-South samples 
at the faint-end, we chose the completeness limit of the GOODS-South as the magnitude cut-off for all samples. The solid circles show the 
best-fit values of M», with the errorbars of the K determinations (similar or smaller errors arc derived in I and B). 
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Fig. 2. — The four plots at the top show the SFR as a function of redshift derived from the 1500 A luminosity densities computed from 
the B-selected (blue), I (green) and I+B-selected (black) FDF, and K-selected (red) GOODS-South field. The points are connected by the 
thick lines for clarity. These SFRs are based on a faint-end slope of the LF of -1.07 as derived from the FDF and GOODS data. The dotted 
lines show the effect of assuming a slope of -1.6. The grey-shaded region shows the eff'ect of dust corrections with correction factors between 
5 and 9, followi ng Adclbcrecr & Stcidcl (2000). The grey symbols show the results (taken from the table of Somcrvillc ct al. 2001) of 
IPascarelle et all 11998 , circled crosses), Stcidcl ct al. (1999, open circles), Madau ot al. (1996, open triangles), Madau ct al. 1 1998 , open 
squares), and (taken from Bunker et al. 2004) Iwata ct al. (2003, filled triangles), Giavalisco ct al. I2004a, filled squares), Bouwe ns et alj 
12003b, filled c ircles), Bouwens et al. (2004, hexagonal crosses), Fontana et al. (2003, filled pentagons). Bunker ct al. 12004, open star), 
IBouwens et a l. 12003a, inverted filled triangles). The plots at the bottom show the SFRs of the four catalogues together (left) and the 
SFRs derived considering the contributions of the galaxies brighter than only (right) . 



